Efficient Fermionic One— Loop RG for the 2D Hubbard Model at Van Hove Filling 
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We propose a novel parametrization of the four-point vertex function in the one-loop one-particle 
irreducible renormalization group (RG) scheme for fermions. It is based on a decomposition of the 
effective two-fermion interaction into fermion bilinears that interact via exchange bosons. Besides 
being more efficient than previous A^-patch schemes, this parametrization also reduces the ambiguity 
of introducing boson fields. We apply this parametrization to the two-dimensional (t, t^)-Hubbard 
model using a novel 0-frequency regularization. 



Introduction. In the past decade the one-loop Wilso- 
nian renormalization group (RG) has been extensively 
used to study the weak coupling instabilities of the two- 
dimensional (MO-Hubbard model 0, i, i, i, i, i, 0, i, 
Isl, El, [H, • Successes include the explanation of d- 
wave pairing tendencies from a repulsive interaction. The 
interplay of antiferromagnetic and d-wave superconduct- 
ing instabilities found for small next to nearest neighbor 
hopping — f resembles the behavior of the high tempera- 
ture cuprates near half filling. Avoiding the artificial sup- 
pression of small momentum particle-hole fluctuations by 
regularizing with temperature, the leading instability for 
strong hopping —t' was found to be ferromagnetism at 
Van Hove filling and triplet superconductivity away from 
Van Hove filling [9, 14]. 

The one-loop truncation in the one-particle irreducible 
RG scheme results in an integro-differential equation 
for the four-point vertex function and the self-energy 
The studies mentioned in the last paragraph ne- 
glect the self-energy and the frequency dependence of 
the four-point function. The remaining three indepen- 
dent momenta of the four-point function (one is deter- 
mined by momentum conservation) are numerically dis- 
cretized with N patches in momentum space. In most 
studies these have been chosen along the Fermi surface 
by power counting arguments. Then a system of ordi- 
nary differential equations of size ~ has to be solved 
numerically. Generically for the two-dimensional Hub- 
bard model at low temperature, some coupling constants, 
which describe the vertex function for certain momentum 
combinations, grow large. In order to determine the na- 
ture of the corresponding instability, susceptibilities are 
calculated in the flow by introducing external boson fields 
coupled to appropriate fermion bilinears [15]. 

In this work we develop a novel parametrization of the 
four-point function in the one-particle irreducible RG 
scheme. Guided by the singular momentum structure 
of the right hand side of the flow equation, we identify 
three channels with distinct singular momentum struc- 
tures. In these channels the effective two-fermion inter- 
action is expanded in fermion bilinears that interact via 
exchange bosons. These bosons, however, are dealt with 
in a purely fermionic language. We show that only a 
small number of terms is needed to capture the essen- 



tial features of the one-loop RG flow. This improves 
the efficiency of the parametrization compared to previ- 
ous A/'-patch schemes. A similar reduction has rece ntly 
been used for the single impurity Anderson model [l6[. 
We apply the method to the two-dimensional Hubbard 
model at Van Hove filling and compare the results with 
the temperature RG flow [9] . Instead of regularizing with 
temperature, we introduce a novel frequency regulariza- 
tion scheme with scale parameter Q. 

Decomposition of the Effective Interaction. Although 
the initial vertex function of the Hubbard model at high 
scales is a constant in momentum space, according to the 
RG equation a non-trivial momentum structure evolves 
in the flow as the scale is lowered. The instabilities of 
the RG flow, are however mainly determined by the sin- 
gular momentum structure of the right hand side of the 
RG equation. Therefore, we develop a parametrization 
of the vertex function that simplifies the momentum de- 
pendence but keeps track of all possible singular contri- 
butions. As long as the vertex function is still regular, the 
only momentum dependence that can change the singular 
behavior of the RG equation is the transfer momentum 
that flows through the scale-derivative of the particle- 
particle or particle-hole two-fermion bubble. 

On the right hand side of the RG equation three dif- 
ferent classes of graphs with distinct transfer momenta 
contribute - the particle-particle graph, the crossed 
particle-hole graph, and the direct particle-hole graphs 
|15| . Accordingly, we identify three different channels, 
which we first assume as general charge and spin rotation 
invariant two-fermion interactions. They are defined by 
three conditions: (1) each vertex function of the chan- 
nels absorbs one singular momentum, that is, one trans- 
fer momentum; (2) the channels describe the interaction 
of Cooper pairs, spin operators, and density operators, 
respectively, if the transfer momentum becomes singu- 
lar; and (3) each channel separately satisfies particle- 
hole symmetry and the fermionic antisymmetry. Then 
the evolution of the superconducting channel is given by 
the particle-particle graph, of the magnetic channel by 
the crossed particle-hole graph, and of the forward scat- 
tering channel by a combination of crossed and direct 
particle-hole graphs, see Figure [TJ The vertices in Fig- 
ure [T] depict the fermionic irreducible four-point vertex 



2 



/c2, /cs), where we denote momentum k and fre- 
quency k() together by k = (k, fep). Spin is conserved 
along the fermion hues, see 



^sc{ki,k3,ki+k2) = 



point symmetry group of the lattice. Then, the expan- 
sion involves non-diagonal terms 




FIG. 1: Evolution of the superconducting, magnetic, and for- 
ward scattering channel respectively. The dot stands for a 
derivative with respect to the scale Q and the doublesided ar- 
row in the last graph symbols the exchange of two momenta 
and frequencies. 

The initial condition for the channels is zero, so they 
are generated in the flow. The full vertex function of the 
Hubbard model is then given by the sum 

Vn{ki,k2,k3) = U -^%{ki,k3,ki^k2) 

+ ^^(fci, fe, - fci) + fe, k2 - ks) (1) 

- ^^^{kuk2,k2-k3). 

This decomposition is exact on one-loop level since all 
graphs with their full momentum and frequency depen- 
dence are taken into account. The assignment of the 
graphs, however, can be ambiguous for special momenta. 
Using the Fierz identity, for example, a constant term 
can be freely distributed among the channels. This is 
why the initial Hubbard repulsion U is not decomposed 
and kept explicitly in the decomposition as a constant. 
However, this does not mean that the on-site term in 
the full interaction remains independent of scale. The 
corrections to it are rather absorbed as contributions to 
the above defined three channels. 

Next we further decompose the two-fermion interac- 
tion by writing each channel as a sum of terms, in 
which two fermion bilinears interact via an exchange 
boson propagator, where the boson momentum is the 
corresponding singular transfer momentum of the chan- 
nel. For fixed transfer momentum / the vertex function 
^^^{q, q^l) = ^^^{q\ q^l) can be seen as the kernel of 
a Fredholm operator. Therefore, due to compactness, a 
diagonal expansion in eigenfunctions is possible in prin- 
ciple. To avoid a calculation of the (scale-dependent) 
eigenfunctions, we simply expand in an ort honor mal ba- 
sis of TT, tt]^) of form factors fn that represent the 
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RU<l,q',l). 



(2) 



For simplicity the form functions are chosen frequency 
independent here. The expansion coefficients D^^(l) de- 
pend on the transfer momentum and frequency / and are 
called superconducting boson propagators. In practise, 
we restrict to a finite set X of terms. The remainder 
function R^^ accounts for the error. 

Similarly, the magnetic and forward scattering chan- 
nels are expanded in appropriate frequency and scale in- 
dependent form factors 

*m(9,'/',0= E ^m„(0/m(q+^)/«(q'-^) 

+ Rl{q,q',l) (3) 

^^{q,q\l)= K^n{l)fm{ci+\)fn{c{-\) 

+ R^{q,q',l). (4) 

This introduces the magnetic boson propagators M^^, 
the forward scattering boson propagators i^^^, and ad- 
ditonal remainder terms. 

The Flow of Boson Propagators. The fiow equations 
for the boson propagators are obtained by inserting the 
proposed decomposition in the RG equation for the ver- 
tex function and projecting each channel onto the 
coefficients of the orthogonal form factor expansion. If 
the remainder terms are dropped, this gives a closed sys- 
tem of integro-differential equations for the boson prop- 
agators. Compared to the original RG equation for the 
vertex function, their solution is now given by several 
functions dependent on one momentum and frequency 
instead of one function dependent on three momenta and 
frequencies. This is a considerable simplification and, in 
particular, makes the numerical implementation more fa- 
vorable. Furthermore, like the two-fermion bubbles, the 
boson propagators have only point singularities in mo- 
mentum space instead of extended Fermi surface singu- 
larities. 

In the case of a regular Fermi surface with irrelevant 
Umklapp scattering the RG fiow is dominated by the su- 
perconducting boson propagators if the initial interaction 
contains an attractive part. This scenario is well under- 
stood even analytically and there are clear arguments 
which form factors have to be chosen, such that the re- 
mainder terms are negligible. If the initial interaction is 
purely repulsive, we are able to show that contributions 
of the magnetic channel induce an attractive d-waye pair- 
ing interaction in the superconducting channel. 

At Van Hove filling, however, the Fermi surface has 
saddle points, which cause a logarithmic singularity in 
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the density of states. Furthermore, near half fihing, Umk- 
lapp scattering processes become relevant. Then there is 
strong mixing between particle-particle and particle-hole 
channels. This case is not analytically understood yet 
and treated numerically in the following. Here the form 
factors are chosen guided by previous A/'-patch studies 
and the remainder terms are neglected. A rigorous anal- 
ysis of the remainder terms is left for future work. 

The numerical implementation of the boson propaga- 
tor flow is eased by specifying the momentum and fre- 
quency dependence of the boson propagators. Here we 
first neglect the frequency dependence of the boson prop- 
agators. Together with the frequency independent form 
factors, this corresponds to a vertex function that does 
not depend on frequency at all. The right hand side of 
the flow equations is evaluated at boson frequency zero, 
which gives the main contribution. An analysis of the 
two-fermion bubbles indicates that the boson propaga- 
tors can only become singular at boson momentum (0, 0) 
and (7r,7r). Therefore we split up each boson propaga- 
tor into a part around (0,0) and a part around (7r,7r). 
In these regions we parameterize the momentum depen- 
dence of the boson propagators by step functions, postu- 
lating a higher accuracy at the singular points. For ex- 
ample, a peak of the magnetic boson propagator Mss(l) 
with constant form factor /s(q) = 1 at 1 = (0,0) indi- 
cates ferromagnetic and at 1 = (tt, tt) antiferromagnetic 
fluctuations. 

In order to set up the RG flow, a regularization pa- 
rameterized by an RG scale has to be introduced. Since 
we are especially interested in the interplay between d- 
wave superconductivity and ferromagnetism, we have to 
choose a regularization that is sensitive to small momen- 
tum particle-hole fluctuations [9]. Here we introduce a 
frequency regularization by multiplying the bare propa- 
gator with the function 

(5) 

where po is the frequency part of p and Q is the RG scale 
frequency. For 11 > all graphs are regularized in the 
infrared and for Q ^ the original model is recovered. 
Unlike a Fermi surface cut-off, this regularization repro- 
duces the correct (In (1)^ -scaling of the particle-particle 
and In ll-scaling for the particle-hole bubble at Van Hove 
filling. Compared to the temperature flow, this regular- 
ization allows a clear definition of the initial condition of 
the RG flow. We start the flow at a high initial scale 
and treat the scales 1] > l^o by perturbation theory in 
the coupling constant /7 > to second order. Due to the 
infrared regularization, perturbation theory converges if 
U/Qq is small enough. For large enough Qq^ we find that 
the dependence on is negligible. 

Results and Outlook. Here we present results ob- 
tained by expanding each channel with two form fac- 
tors, namely isotropic 5-wave /s(q) = 1 and dx2_y2- 
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FIG. 2: The critical scale Qc in dependence on next to nearest 
neighbor hopping —t' at temperature zero with initial scale 
Qo = 15t. The chemical potential is set to match Van Hove 
filling at each — 1^ The instabilities of the Landau Fermi 
liquid are determined as antiferromagnetism (AF), dx2_y2- 
superconductivity (dSC), and ferromagnetism (FM). 



wave /i(q) = cosgi — cosg'2. Amongst others, the corre- 
sponding boson propagators allow to detect instabilities 
towards s- and dx2_y2-wdiYe superconductivity, ferro- 
and antiferromagnetism, forward and exchange scatter- 
ing, and also a Pomeranchuk instability. It turned out, 
however, that the boson propagator flow applied to the 
(t, t')-Hubbard model at zero temperature and Van Hove 
filling already gives reasonable results if only the constant 
form factor is taken into account in every channel plus the 
dx'2-y2-wdive form factor in the superconducting channel. 
Although 5-wave superconductivity is suppressed by the 
initial repulsion, its screening effect is essential for the 
RG flow. 

Starting from U = 3t we observe a generic flow to 
strong coupling. If the maximum of one boson propaga- 
tor reaches the value 20t, the flow is manually stopped at 
the thus defined "critical" scale Qc Like in the previous 
A/'-patch studies of the symmetric phase, this is inter- 
preted as an instability towards a corresponding ordered 
state. The critical scale is plotted over next to nearest 
neighbor hopping —f/t in Figure [21 where from now on 
t = 1. For small —t\ where Van Hove filling is close to 
half filling, antiferromagnetism dominates. Since the flow 
is stopped at a relatively high scale, perfect nesting-like 
effects are not restricted to —f = 0. If —f is increased, 
we observe a tendency to incommensurate antiferromag- 
netic order (marked with open squares). For intermediate 
—t' the leading instability is 6/3.2 _ ^2 -superconductivity. 
It is induced by antiferromagnetic fluctuations, as can 
be seen from the flow equations. For high —t' ferro- 
magnetism is the dominant instability. In the interval 
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—t' G (0,0.34) the critical scale drops by two orders of 
magnitude. 

The results obtained with the decomposed RG flow 
agree well with the temperature TV-patch flow [9]. The 
proposed approximation to the one-loop RG flow re- 
produces the qualitative instabilities and the values for 
—t', where the transitions between different instabilities 
take place. Indicated by the decreased critical scale, 
(i3,2_ ^2 -superconductivity and ferromagnetism suppress 
each other. For comparison, the Stoner criterion for 
ferromagnetism, which is obtained by neglecting the 
particle-particle channel, is plotted in Figure[2]on the fer- 
romagnetic side. There is one significant deviation from 
the temperature flow: in our present scheme, the sup- 
pression of Qc in the transition region from d-wave su- 
perconductivity to ferromagnetism is much weaker. Since 
both the temperature flow and our calculation involve ap- 
proxiamtions, the existence of a quantum critical point 
remains an open question. 

In summary, the proposed decomposition of the effec- 
tive two-fermion interaction is efficient for studying com- 
peting instabilities in the (t, t')-Hubbard model. Sep- 
arating leading from subleading processes reduces the 
complexity of the one-loop flow equations. Although 
we cannot yet rigorously justify dropping the remainder 
terms at Van Hove filling, there are clear arguments for 
this for regular Fermi surfaces. In fact, the comparison 
with previous A/'-patch studies shows that the qualitative 
structure of the RG flow is preserved. 

In the numerical implementation the momentum de- 
pendence of the boson propagators is discretized using 
step functions. This is more precise than the general 
patching in A^-patch schemes since the choice of step 
functions can be guided by the one-loop bubbles. Gen- 
erally, if the momentum dependence of the one-loop bub- 
bles can be is parameterized in an analytical form, then 
it is possible to extract a functional parametrization of 
the boson propagators from the flow equations, at least 
for small momenta. Similarly, the dependence on small 
frequencies can be taken into account. Deviations from 
the large frequency behavior would then be subject to 
further remainder terms. 

The decomposition of the effective two-fermion inter- 
action into sums of fermion bilinears interacting via bo- 
son propagators allows to decouple the fermion bilin- 
ears by multiple Hubbard Stratonovich transformations. 
Thereby the ambiguity of introducing boson fields is not 



completely removed as discussed above. However, due to 
the definition of the channels based on the singular mo- 
mentum structure, it is reduced. Thus our results serve 
as an improved starting point for a continuation of the 
RG flow into the symmetry broken phase in a (partially) 
bosonized form [ItI, [H, [11, 0, EH . 
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